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ABSTRACT 

We present 0.5-2 keV, 2-8 keV, 4-8 keV, and 0.5-8 keV (hereafter soft, hard, ultra-hard, and full bands, respectively) 
cumulative and differential number-count (log N -log S ) measurements for the recently completed ~4 Ms Chandra 
Deep Field-South (CDF-S) survey, the deepest X-ray survey to date. We implement a new Bayesian approach, 
which allows reliable calculation of number counts down to flux limits that are factors of rH .9-4.3 times fainter 
than the previously deepest number-count investigations. In the soft band (SB), the most sensitive bandpass in 
our analysis, the ~4 Ms CDF-S reaches a maximum source density of ~27,800 deg 2 * . By virtue of the exquisite 
X-ray and multiwavelength data available in the CDF-S, we are able to measure the number counts from a variety 
of source populations (active galactic nuclei (AGNs), normal galaxies, and Galactic stars) and subpopulations 
(as a function of redshift, AGN absorption, luminosity, and galaxy morphology) and test models that describe their 
evolution. We find that AGNs still dominate the X-ray number counts down to the faintest flux levels for all bands 
and reach a limiting SB source density of R» 14,900 deg -2 , the highest reliable AGN source density measured at 
any wavelength. We find that the normal-galaxy counts rise rapidly near the flux limits and, at the limiting SB 
flux, reach source densities of ~ 1 2,700 deg 2 and make up 46% ± 5% of the total number counts. The rapid rise 
of the galaxy counts toward faint fluxes, as well as significant normal-galaxy contributions to the overall number 
counts, indicates that normal galaxies will overtake AGNs just below the ~4 Ms SB flux limit and will provide a 
numerically significant new X-ray source population in future surveys that reach below the «s4 Ms sensitivity limit. 
We show that a future ~ 1 0 Ms CDF-S would allow for a significant increase in X-ray-detected sources, with many 
of the new sources being cosmologically distant (z > 0.6) normal galaxies. 

Key words: cosmology: observations - galaxies: active - galaxies: starburst - X-rays: galaxies 
Online-only material: color figures, Supplemental data file (tar.gz) 


1. INTRODUCTION 

Deep extragalactic X-ray surveys conducted with Chandra 

and XMM-Newton have resolved the vast majority of the 

0.5-10 keV cosmic X-ray background (CXRB) and have pro- 

vided substantial new insight into the X-ray point sources de- 
tected (see, e.g., Brandt & Hasinger 2005; Brandt & Alexander 

2010). A fundamental quantity used to characterize the extra- 

galactic X-ray source population is the cumulative X-ray num- 
ber counts, which quantify how the cumulative number of X-ray 

sources per unit area, N, increases with decreasing flux, S (also 
referred to as log fV-log S; see, e.g., Brandt et al. 2001; Rosati 

et al. 2002; Moretti et al. 2003; Bauer et al. 2004; Kim et al. 

2007; Georgakakis et al. 2008). The number counts provide 

an observational constraint that must be considered when con- 

structing physical models that describe the X-ray evolution of 

extragalactic sources in the universe. For example, successful 

models including the supermassive black hole (SMBH) accre- 
tion history of the universe (e.g., Gilli et al. 2007; Treister et al. 

2009) and the X-ray evolution of normal galaxies (primarily 

driven by X-ray binaries and hot gas; see, e.g., Ranalli et al. 


2003, 2005) must predict X-ray source densities consistent with 
the observed number counts. 

Due to its high angular resolution and low background, 
Chandra is the only current X-ray observatory capable of 
providing new constraining data in the ultra-deep regime 
(below 0.5-8 keV fluxes of a few x 10~ 16 erg cm -2 s -1 ), since 
the deepest Chandra surveys have now greatly surpassed the 
XMM-Newton confusion limit (see, e.g., Brandt et al. 2001; 
Giacconi et al. 2001, 2002; Alexander et al. 2003; Luo et al. 
2008; Xue et al. 2011). As new Chandra surveys continue to 
probe the extragalactic X-ray universe to fainter depths, the 
number counts continuously rise as fainter X-ray populations 
are revealed. Naturally, it is at the faintest flux levels where 
ultra-deep Chandra surveys are probing new regions of discov- 
ery space and classes of extragalactic sources that were poorly 
sampled at X-ray energies in the past (e.g., star-forming and 
passive galaxies, and obscured and low-luminosity active 
galactic nuclei (AGNs)). 

Bauer et al. (2004; hereafter B04) measured the number 
counts for X-ray point sources detected in the %2 Ms Chandra 
Deep Field-North (CDF-N) and «1 Ms CDF-South (CDF-S), 
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the deepest Chandra surveys at the time. Using the avail- 
able multiwavelength data, B04 were able to distinguish 
between different X-ray-emitting populations (e.g., AGNs, 
normal galaxies, and Galactic stars) and measure their con- 
tributions to the CXRB. B04 found that, generally, AGNs make 
up exclusively the relatively bright number counts. However, 
at the faintest flux levels in the 0.5-2 keV bandpass, the most 
sensitive bandpass studied in the CDF surveys, B04 showed that 
normal galaxies start to compose an appreciable fraction of the 
number counts («25% of sources with So.5-2keV ^ (2-10) x 
10“ 17 erg cm" 2 s _1 ). 

Recently, the deepest extragalactic X-ray survey yet con- 
ducted, the CDF-S, has reached a total exposure of ~4 Ms (Xue 
et al. 2011, hereafter XI 1). XI 1 presented point-source catalogs 
and data products for the survey and provided basic multiwave- 
length classifications for the 740 individually detected X-ray 
sources in their main catalog. This work has revealed that, at 
the faintest flux levels, normal galaxies are playing an increas- 
ingly important role in the new sources detected. For example, at 
0.5-2 keV fluxes below «5 x 10“ 17 erg cm" 2 s _1 , »50% of the 
X-ray-detected sources are classified as likely normal galaxies; 
below the current detection limits, it is almost certain that the 
normal-galaxy fraction continues to increase (e.g., B04; Ranalli 
et al. 2005). 

In this paper, we present number counts for the new ~4 Ms 
CDF-S, focusing on the faint-flux regime (0.5-2 keV fluxes 
<10“ 17 erg cm" 2 s“*). In Section 2, we describe a new Bayesian 
method for computing number counts, which properly accounts 
for biases and measurement uncertainties that are present in the 
important flux regime near the detection limit. In Section 3, we 
utilize the available multiwavelength photometric data as well 
as optical spectroscopic and photometric redshift catalogs to 
quantify, with good reliability, the relative contributions to the 
number counts from AGNs, normal galaxies, and Galactic stars. 
We further break down the number counts to quantify the con- 
tributions that subpopulations make to the AGN (with redshift, 
intrinsic AGN column density, and X-ray luminosity subpop- 
ulations) and normal galaxy (with redshift and morphological 
subpopulations) number counts. We conclude Section 3 by dis- 
cussing how each population and subpopulation contributes to 
the overall CXRB intensity. In Section 4, we use our subpopula- 
tion number-count estimates to compare with predictions from 
phenomenological models that describe how the X-ray emis- 
sion from accreting SMBHs and normal galaxies are expected 
to evolve with cosmic time. In Section 5, we use these models 
to estimate directly number counts to flux levels below current 
detection limits and highlight the prospects of deeper X-ray 
surveys. In Section 6, we summarize our findings. 

Throughout this paper, we make use of the main point-source 
catalog and data products provided by X 1 1 . The Galactic column 
density for the CDF-S is 8.8 x 10 19 cm" 2 (Stark et al. 1992). 
All of the X-ray fluxes and luminosities quoted throughout this 
paper have been corrected for Galactic absorption. In the X-ray 
band, we make use of four bandpasses: 0.5-2 keV (soft band 
(SB)), 2-8 keV (hard band (HB)), 4-8 keV (ultra-hard band 
(UHB)), and 0.5-8 keV (full band (FB)). Values of //„ = 
70 km s _1 Mpc“*, £2 M = 0.3, and £2 A = 0.7 are adopted 
throughout this paper (e.g., Spergel et al. 2003). 


flux levels. In the faint-flux regime, computing number counts 
presents a challenge, since these calculations must properly ac- 
count for (1) the non-negligible spatial variations in sensitivity 
across the Chandra image, (2) incompleteness issues related to 
source-detection algorithms (see Section 2.1), and (3) the Ed- 
dington bias. The B04 investigation of the CDF number counts 
implemented Monte Carlo simulations using an extrapolated 
faint-flux number-count model to measure and correct for these 
biases. Although this method provides reasonable first-order 
corrections that account for the biases near the flux limits, it does 
not optimize the input faint-flux model. Such an optimization 
would require several Monte Carlo runs with varying faint-flux 
extrapolations, which is not feasible due to the computational 
requirements of this procedure. Such a limitation will therefore 
introduce small systematic errors in the number-count measure- 
ments near the flux limit if the wrong faint-end number-count 
prior is chosen. 

To mitigate this difficulty, we employ an approach similar 
to the Poissonian-based methods described by Georgakakis 
et al. (2008, hereafter G08); however, our method has been 
adapted to account for the different source-detection methods 
adopted for the ~4 Ms CDF-S catalog by XI 1 . The XI 1 source- 
detection criteria make use of both wavdetect (Freeman et al. 
2002) for initial source selection and ACIS Extract (AE) to 
improve photometry and re-evaluate source-detection signifi- 
cance (Broos et al. 2010). AE makes use of binomial statistics to 
evaluate source significance, and it properly accounts for point- 
spread function (PSF) variations between observations and un- 
certainties in local background measurements (see below). Our 
method for computing number counts, described in detail be- 
low, uses a Bayesian approach with maximum-likelihood opti- 
mizations to account for the Eddington bias and completeness 
limitations without requiring a large number of time-consuming 
simulations. 


2.1. Source Recovery Functions and Flux 
Probability Distributions 

Number-count computations at the flux limits depend sensi- 
tively on the solid angle of the survey over which a source of a 
given flux could be detected. As described in Section 3 of XI 1, 
to be included in the main catalog, an X-ray source must (1) be 
detected by wavdetect at a false-positive probability threshold 
of 10“ 5 and (2) contain s counts (derived from AE; see Columns 
8-16 of Table 3 in XI 1) within an aperture representing ~90%' 
of the point-source encircled-energy fraction (EEF) 14 that 
satisfies the following binomial probability criterion: 


p(x > s) = 


x!(n — x)! 


p x ( i - pT 


c P, 


thresh > 


(i) 


where n = s + b ext and p = 1/(1 + b ext /b src ). Here b ext is 
the total number of background counts extracted from a large 
region outside of the point source (while masking out regions 
from other X-ray-detected sources) that was used to obtain 
an estimate of the local background count rate. The quantity 
b SIC is the estimated number of background counts within the 
source extraction region, which was measured by rescaling 


2. METHODOLOGY FOR COMPUTING 
NUMBER COUNTS 

Our primary goal is to evaluate the number counts across 
the entire »M65 arcmin 2 ~4 Ms CDF-S to the faintest possible 


14 As noted by XI 1, in special cases where sources had nearby neighbors, the 
source-detection probability was measured using an aperture smaller than the 
^90% EEF; however, only ss2%— 3% of all sources had moderate source 
crowding, where detection probabilities were measured using apertures that 
encompass <50% EEF. We therefore expect that this will have a negligible 
effect on our number-count measurements. 
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b exl to the area of the source aperture (typically, b ext /b SIC ~ 
16; see Section 4.1 of XI 1 for details). In Equation (1), 
^thresh = 0.004, the value adopted by XI 1. This choice of 
P thresh was empirically selected both to optimize the number 
of sources detected and to ensure that nearly all detected 
sources are reliable (see Section 4.1 of XI 1). This multi- 
stage procedure for identifying a highly reliable list of source 
candidates will not result in a complete selection of all real 
sources in the image having binomial probabilities ^ /Wsh, 
since the initial wavdetect selection has more complex source- 
detection criteria (see Freeman et al. 2002 for details) than 
the simple criterion given in Equation (1) (see Section 2.2 and 
the Appendix for discussion on correcting for these issues when 
computing number counts). 

Using a probabilistic approach, we compute the solid angle 
of the sky within which a source with intrinsic flux S would 
be detected if present. First, using Equation (1), we define 
the quantity L, which is the number of counts required for a 
detection, which satisfies the relation P(x Js L) = 

^thresh — 

0.004. In order to compute L at each location in the CDF-S 
image, it is necessary to estimate local values of b src and b ext that 
are similar to those that would be estimated following the XI 1 
approach used for the X-ray-detected sources. Unfortunately, the 
X 1 1 procedure for extracting these values at every location in the 
image is computationally prohibitive, since source photometry 
is performed on an observation-by-observation basis (i.e., in 
up to 54 observations in total) and makes use of polygonal 
apertures that approximate local PSFs (see XI 1 for further 
details). To overcome this issue, we estimated b SIC by extracting 
background counts from the merged CDF-S background maps 
(see Section 7.1 of XI 1) using circular apertures with sizes that 
encompass the 90% EEF for a point source. This approach was 
tested by comparing our b STC values measured in the regions of 
the XI 1 sources with those of XI 1. We find good agreement 
between values and a la scatter at the «U8%-25% level. To 
estimate appropriate values of b ext at each location (pixel) on 
the CDF-S image, we measured the local off-axis angle 6 p ; 
using the XI 1 point-source catalog, we adopted the maximum 
value of b ext for sources with off-axis angles 0 = 0 P ± 0(25. 
Given values of b SIC and b ext , we numerically solved the relation 
P( x ^ L) = 0.004 to obtain L at each image location. In this 
manner, we constructed a spatial sensitivity map consisting of 
L values. 

In principle, L could be combined with the exposure time to 
estimate a count-rate limit, which can in turn be converted to a 
single flux limit using a count-rate-to-flux conversion (as done 
in the XI 1 catalogs); however, this single flux-limit approach 
does not directly incorporate the probabilistic nature of source 
detection that is important in the low-count regime. It is therefore 
more informative to compute the probability that a source with 
intrinsic flux, S, would be detected given that L counts are 
required. Such a source is expected to contribute the following 
number of counts to the source-detection cell: 

T = t exp Cr/S + b src , (2) 

where t ex p , C, and i] are the effective exposure time, the 
conversion from flux to count rate, and the EEF, respectively. 
Values of f cxp are taken directly from the exposure maps from 
XI 1 (see their Section 3.1). These exposure maps include the 
effects of vignetting, gaps between CCDs, bad-column filtering, 
bad-pixel filtering, and spatial and time-dependent degradation 
in quantum efficiency due to contamination on the ACIS optical- 
blocking filters. Values of C will depend on the spectral energy 



S = 0.5-2 keV Flux (ergs cm 2 s 

Figure 1 . Detection probability as a function of intrinsic 0.5-2 keV flux for 
point sources located on axis (dashed curve) and ~8' off axis (solid curve). 
These curves were computed considering a counts-to-flux conversion factor 
appropriate for a power-law SED with T = 1 .4 (see Equations (2) and (3)), 
the mean photon index of the CXRB. Each curve asymptotes to the defined 
detection probability threshold Pthresh = 0.004 (dotted line) at low flux levels. 
To construct the cumulative solid angle as a function of flux limit for the entire 
ss4 Ms CDF-S image, we computed such curves at each location on the Chandra 
image multiplied by the effective solid angle per pixel (RiO.242 arcsec 2 ) and 
added them (see Figure 2). 


distribution (SED) shape. In this work we use power-law 
SEDs to characterize the observed spectra; these SEDs can 
be described using the photon index U For each source, T 
was derived using the HB-to-SB count-rate ratio (corrected 
for differential vignetting and exposure times) as a proxy for 
spectral slope. The HB-to-SB count-rate ratio was calibrated 
against T using the AE-automated XSPEC-fitting procedure for 
relatively bright X-ray sources (with FB counts greater than 200; 
this ensures reliable XSPEC-fitting results). This approach takes 
into account the multi-epoch Chandra calibration information 
(see XI 1 for further details). 

Therefore, at each location on the image, the probability of 
source detection as a function of flux S can be computed as 
follows: 


L+bpy 


Pdet — ^2 


( L + b ex t )! 


x=L 

x I 1 


x\{L + b ext - x)l V T + b, 

L+b ca -x 


T + b„ 


( 3 ) 


where b ext is the extracted background counts that were used 
to calculate the local background in Equation (1). Figure 1 
shows P det as a function of 0.5-2 keV flux -So. 5-2 keV for 
randomly selected positions at 9 ss O' and 8 ' offsets from 
the average aim point of the 4 Ms CDF-S (assuming ri — 0.9 
and a counts-to-flux conversion factor appropriate for a power- 
law SED with r = 1.4). We note that as the flux drops to 
zero, P d et asymptotically approaches our detection threshold 
Pthresh = 0.004. This shows that, even when no source is present, 
there is still a finite (yet small, i.e., P = 0.004) probability that a 
positive fluctuation may exceed our adopted detection threshold; 
however, as shown by X 1 1 through multi wavelength counterpart 
matching, the initial wavdetect source selection ensures that 
very few false sources are present in the XI 1 catalog. 

We consider that each Chandra pixel represents a small local 
solid angle of size dCl & 0.242 arcsec 2 (i.e., «0'.'492 x O'.' 492) 
with a flux and SED (i.e., E) dependent detection probability 
distribution characterized by Pjet- Under this consideration, the 
total surveyed solid angle over which sources with flux S and 
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Figure 2. Effective solid angle A(S) vs. intrinsic flux limit S for the SB 
(dotted curve), HB (short-dashed curve), UHB (dash-dotted curve), and FB 
(solid curve). These curves have been computed for a power-law SED with 
r = 1.0, the median photon index for sources detected in both the SB and 
HB. The shaded regions give the effective solid angle curves appropriate for 
the interquartile (&25%-75%) range F = 0.4-1 .5. These curves were computed 
following the probabilistic methods discussed in Section 2.1 and account for 
uncertainties in measured flux conversions at the lowest flux levels where the 
number of detected counts is small. The horizontal dotted line represents the 
«JlO arcmin 2 limit, above which we can confidently compute number counts 
(see Section 2.1 for details). These curves are qualitatively different in nature 
from those produced by XU, which utilize single-valued count-rate-to-flux 
conversions. 


photon index T could be detected when present is therefore 
A(S, F) = / J dct.; dU. (hereafter the effective solid angle), 
where the summation is over all possible detection cells (i.e., 
all pixels). In Figure 2, we show A(S, T = 1.0) versus S 
for the «s4 Ms CDF-S in the four bandpasses. The shaded 
regions show A(S, T) in the range of V — 0.4-1. 5 (median value 
r = 1.0), which represents the interquartile (i.e., 25%-75%) 
range for the 332 XI 1 sources with estimates of F that were 
not based on limits. We note that generally A(S, T) increases 
with increasing T. We find that the A(S, T) curves asymptote 
to a value of ~2 arcmin 2 approaching .S' = 0, suggesting 
that sources with extremely low fluxes (down to zero) could 
in principle produce fluctuations exceeding the probability 
threshold defined in Equation (1). However, as discussed above, 
such sources would most efficiently be removed using our initial 
wavdetect screening. Therefore, we cannot use information 
at such flux levels to determine number counts reliably. We 
therefore choose to restrict our number-count computations to 
flux levels where >10 arcmin 2 solid angle is accessible in our 
survey for the case of T = 1.4 (the mean SED of the X-ray 
background; e.g., Moretti etal. 2009). Our adopted & 1 0 arcmin 2 
solid-angle limit additionally constrains our flux limits. The 
resulting limits are 5.1 x 10~ 18 , 3.7 x 10~ 17 , 4.6 x 10~ 17 , and 
2.4 x 10“ 17 erg cm 2 s" 1 for the SB, HB, UHB, and FB, 
respectively; these are factors of ^1.9-43 times fainter than 
those of B04 and G08, the previously deepest number-count 
studies. For the full range of X-ray spectral slopes (i.e., T values) 
in the XI 1 sample, ~5-20 arcmin 2 of solid angle is available 


for number-count computations at the flux limits (above the 
asymptotic regime). We note that the flux limits derived here are 
fainter than those presented in XI 1, which were ~9. 1 x 10~ 18 
and ^5.5 x 10 -17 erg cm~ 2 s -1 in the SB and HB, respectively 
(see also the area curves in our Figure 2 compared with 
Figure 23 in XI 1). This is due to the fact that XI 1 considered 
only a single count-rate-to-flux conversion factor and did not 
use the probabilistic approach adopted here. 

To account for the fact that, for each X-ray-detected source 
in our main catalog, we are only able to measure reliably 
the total observed counts s (see above), the local background 
b exU and that the intrinsic flux S may be subject to large 
uncertainty (particularly in the low-count regime), we consider 
the conversion from counts to flux for each source to be 
probabilistic. For each X-ray-detected source, we computed the 
flux probability distribution as follows: 


(s+b ext y. / T 
j!6 ext ! \T + b ex t 


S 


X 


T + b e 


dN/dS | 


model ’ 


(4) 


where the term dN /dS\ moAA is a Bayesian prior, based on the 
differential number counts, which accounts for the Eddington 
bias near the sensitivity limit. As noted in B04, the slope of the 
number counts of AGNs, normal galaxies, and Galactic stars will 
differ at the flux limit of the CDF-S. Therefore, dN/dS | model 
for a given source will depend on which source population it 
belongs to. Previous studies of X-ray number counts (e.g., Rosati 
et al. 2002; B04; Kim et al. 2007; G08) have shown that power 
laws provide good fits to the overall shapes of the log /V-log S. 
To first order, we use priors based on the following power-law 
parameterizations : 

d_N_ AW _ 
dS 

v JK AGN (S/S ref )-^ GN (S^/ b AGN ) 

U AGN (/ b /S ref )ft AGN -^° N ( 5 / Sref )-ft AGN (5 > /AON) 


^ Sal = K gal (S/S ref )-^ 
dS 

/7 N star 

^ = K^(S/S ret y^, ( 5 ) 

dS 

where dN/dS( AGN), dN/dS( gal), and dN /dS(star:) are differ- 
ential number-count parameterizations to be applied to AGNs, 
normal galaxies, and Galactic stars, respectively, / b AGN is the 
flux related to the break in the double power law used to de- 
scribe the AGN number counts, and S ie f = 1 0~ 14 erg cm~ 2 s _1 . 
As we show in Section 3.1 below, we characterize each 
X-ray source using the X-ray and multiwavelength data and 
provide best estimates of the parameters in Equation (5) (i.e., 
K, ft, and f, values) for the AGNs, normal galaxies, and 
Galactic stars. For each source, we computed Pj down to the 
flux limits defined above and normalize Equation (4) using 
f P T dS = 1. In Figure 3, we show examples of Pj as a 
function of 0.5-2 keV flux for two AGNs in the main Chandra 
catalog having AE probability P ~ 10 47 (solid curve) and 
P ~ 0.003 (dashed curve; near the detection threshold). For this 
computation, we utilized values of ft AGN = 1 .49, ft AGN = 2.49, 
and / b AGN = 5.6 x 10 -15 erg cuT 2 s -1 (the normalization is 
arbitrary); as we will show in Section 3.1, these values represent 
the best-fit parameterizations for AGN SB number counts. 
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S = 0.5-2 keV Flux (ergs cm 2 s ’) 

Figure 3. 0.5-2 keV flux probability distributions for two AGNs in the XI 1 
catalog with no-source-detection probabilities (provided by AE) of ~10~ 45 
(solid curve) and Si0.003 (dashed curve). The latter source is roughly at the 
boundary of our adopted detection threshold Pthresh = 0.004 and therefore 
illustrates the extent to which the flux conversion for sources in the low-count 
regime with non-negligible Eddington biases can be affected. The vertical dotted 
lines indicate the most probable single- value flux given in the XI 1 catalog. 


2.2. Cumulative Number-count Computation 

Using the above information, we computed cumulative num- 
ber counts for the ~4 Ms CDF-S survey using the following 
integral: 


N{> S ) = 


-[ 


Y'PTjWt/AiSSi) 


dS, 


(6) 


_ /= 1 


where the summation is over all X-ray-detected sources in a 
particular category of sources (e.g., all sources, AGNs, normal 
galaxies, etc.) and F, represents the effective photon index for 
source i. For the 332 X-ray sources detected in both the SB and 
FIB, we utilized the photon indices from XI 1; otherwise, we 
used r = 1.4. Since Pjj is dependent on our input Bayesian 
prior (which depends on source type), the number counts will 
depend mildly on the choice of our model. The term C' is a count 
rate and off-axis angle-dependent completeness term, which 
corrects for the fact that our X-ray point-source catalogs were 
constructed by (1) running wavdetect to form an initial list of 
candidate sources and (2) assessing the probability of detection 
using AE. If we chose to use a catalog based strictly on the 
AE source-detection probability assessment for every location 
on the image, then C' = 1 for all sources; however, such a 
catalog would produce large numbers of false sources, leading 
to significant errors in the number counts at faint fluxes, and 
is computationally impractical. In practice, the completeness 
corrections are small perturbations (median C' 1.00 for 
all bands, with «95% of sources having C' < 1. 7-2.4) and 
affect number counts near the flux limits. C' and its count rate 
and off-axis angle dependencies were measured using marx 
simulations, and a full description of this procedure can be 
found in the Appendix. 

Following the same approach implemented in Section 5 of 
G08, we utilized maximum-likelihood techniques to optimize 
the Bayesian model parameters. As discussed above, we made 
use of power-law differential number-count models for our 
Bayesian priors (see Equation (5)), which are dependent on 
source type. For a given source, we first characterized its type 
(i.e., AGNs, normal galaxy, or Galactic star; see Section 3.1 
below), chose an appropriate model in Equation (5), and 
computed the probability p, of that source being present in the 


~4 Ms CDF-S catalog as follows: 

= f P T {S, Nj )d S 
P ‘ fdN/dSl^A&WS’ 

where Nj represents the total number of counts measured within 
the extraction cell for each detected source. Considering our 
model, the total likelihood of obtaining the ~4 Ms source catalog 
and its source-count distribution can be computed as ]~[ ( . . We 

maximized the total likelihood for the model input parameters 
for AGNs, normal galaxies, and Galactic stars separately to find 
best-fit values for each of the parameters in Equation (5). 

3. RESULTS 

3.1. Number Counts by Source Type and Total Counts 

The extensive multiwavelength data in the ?»4 Ms CDF-S 
region allow for the robust characterization of the X-ray- 
detected sources. As described in XI 1, 716 of the 740 
X-ray-detected sources have multiwavelength counterparts, and 
673 have either secure spectroscopic or reliable photometric 
redshifts. Following a similar scheme to that provided by XI 1, 
we have classified the X-ray sources in our sample as AGNs, 
normal star-forming galaxies (hereafter galaxies), and Galactic 
stars. In this process, we make use of six AGN selection criteria. 

1. X-ray luminosity. A source with an observed-frame SB, 
HB, UHB, or FB luminosity (i.e., observed-frame Lx = 
\nd\fx, where dr is the luminosity distance) >3 x 
10 42 erg s -1 is classified as a luminous AGN. In the 
nearby universe, the most powerful star-forming galaxies 
(luminous infrared galaxies) all have Z. 0 . 5-8 keV < 
10 42 erg s _I (e.g., Iwasawa et al. 2009, 201 1; Lehmer et al. 
2010; Pereira-Santaella et al. 2011). 

2. X-ray spectral shape. A source with an effective photon 
index of T < 1 .0 is most likely to be dominated by a single 
powerful X-ray source that is significantly obscured. We 
consider such sources to be obscured AGNs. 

3. X-ray-to-optical flux ratio. A source with an X-ray-to- 
optical flux ratio of log(/ x //«) > 1 (where f x = /o.s- 2 k e v, 
/2-8keV, /4-8keV, or /o.5-8k e v) is considered to have 
X-ray emission significantly elevated compared with 
normal-galaxy processes (e.g., hot gas, low-mass and high- 
mass X-ray binary emission) as traced by stellar emission 
(i.e., /«). Such sources were classified as AGNs. 

4. X-ray-to-radio luminosity ratio. A source with substantial 
excess (i.e., a factor of >5) X-ray emission over the 
level expected from pure star formation (i.e., observed- 
frame SB, HB, FB, or UHB luminosity >5 x (8.9 x 
10 17 Z-i.4ghz); Bauer et al. 2002) is considered to be X-ray 
overluminous compared with the X-ray/star formation rate 
(SFR) correlation, and its X-ray emission is therefore likely 
to be dominated by an AGN. This criterion is similar in 
nature to that of criterion 3 but differs in that it allows for 
a more sensitive classification of a subset of actively star- 
forming galaxies (i.e., detected at 1.4 GHz) using a more 
reliable tracer of the intrinsic galactic SFR that is insensitive 
to dust obscuration. 

5. X-ray variability. Sources found to be variable in the X-ray 
band (with probability >95%) on timescales of months to 
years that also have 0.5-8 keV luminosities greater than 
10 41 erg s' 1 are unlikely to be produced by normal-galaxy 
processes (see Young et al. 2012 for details). These sources 
were classified as AGNs. 


5 


The Astrophysical Journal, 752:46 (23pp), 2012 June 10 


Lehmer et al. 


6. Optical spectroscopy. Finally, sources with optical spectro- 
scopic AGN features such as broad and/or high-excitation 
emission lines are classified as AGNs. These sources were 
identified using the spectroscopic catalogs of Szokoly et al. 
(2004), Mignoli et al. (2005), and Silverman et al. (2010) 
(see XI 1 for details). 

A source is initially classified as an AGN if at least one 
of the above six criteria is satisfied. However, we found that 
five of the sources that were classified as AGNs via the above 
criteria had 1-2 keV and 0.5-1 keV band ratios that were 
consistent with being dominated by hot interstellar gas (see 
Danielson et al. 2012). The initial AGN classification for 
these five sources was based on having high X-ray-to-optical 
flux ratios (criterion 3); however, their X-ray luminosities com- 
pared with their B-band luminosities appear to be consistent 
with hot-gas-dominated galaxies (see Danielson et al. 2012 for 
details). We therefore classified these five objects as normal 
galaxies. In addition to classifying sources directly as AGNs 
and normal galaxies, we searched for direct multiwavelength 
indicators of Galactic stars. Likely Galactic stars were identi- 
fied using (1) optical spectroscopy from Szokoly et al. (2004), 
Mignoli et al. (2005), and Silverman et al. (2010); (2) Hub- 
ble Space Telescope ( HST) stellarity indices >0.7 (from the 
Caldwell et al. 2008; GEMS catalogs); and (3) best-fit stellar 
templates from MUSYC photometric redshift fits (Cardamone 
et al. 2010). All Galactic star candidates were visually screened 
and obvious normal galaxies were noted. The remaining stellar 
candidates were then classified as “Galactic stars.” Remaining 
sources that were not classified as either AGNs or stars were 
classified as normal galaxies. Based on these criteria, we there- 
fore estimate that of the 740 X-ray-detected sources, 561 are 
AGNs, 169 are normal galaxies, and 10 are Galactic stars. 

We compared our normal-galaxy classifications with the 
31 radio and X-ray-selected normal galaxies classified by 
Vattakunnel et al. (2012) and found that 25 sources overlap. 
We find that four of the Vattakunnel et al. (2012) sources 
that we classify as AGNs had high X-ray-to-radio luminosity 
ratios (criterion 4), a criterion that was not implemented by 
Vattakunnel et al. (2012). 

In Figure 4, we plot the observed-frame 0.5-8 keV luminosity 
versus redshift for the AGNs and normal galaxies with redshift 
measurements. Our AGN and normal-galaxy samples span 
redshift ranges of z ^ 0.1-8 and z ~ 0.03-2.6, respectively. 
These redshifts are from both photometric and spectroscopic 
redshift estimates, and all redshift estimates above z = 4.76 are 
based on photometric redshifts. We have chosen to use the most 
probable photometric redshift for each source; however, some 
of the z > 4.76 sources have photometric redshifts consistent 
with being at lower redshifts (see Luo et al. 2010 for details). 

Using the classifications adopted above and our methods 
for computing and optimizing number counts described in 
Section 2, we computed the cumulative number counts for 
AGNs, normal galaxies, and Galactic stars and determined the 
best-fit parameters related to their priors (i.e., those character- 
ized in Equation (5)). In Figure 5, we present the breakdown of 
the cumulative number counts from AGNs, normal galaxies, and 
Galactic stars for the four bandpasses (discrete symbols). Total 
number counts for the survey are simply the sum of the number 
counts from each of these populations. In Figure 6, we show the 
differential number counts ( dN/dS ) for the four bandpasses. 
The best-fit dN/dS model parameters, used in computing our 
number counts, have been tabulated in Table 1, and the cumu- 
lative and differential number counts corresponding to these 
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Figure 4. Observed-frame 0.5-8 keV luminosity vs. redshift for AGNs (blue 
triangles) and normal galaxies (red squares). Luminosities were calculated using 
the best redshift estimates and 0.5-8 keV fluxes provided by XI 1 and our 
adopted cosmology. The dashed curve indicates the observed-frame luminosity 
corresponding to the flux limit at the center of the «4 Ms CDF-S. The normal- 
galaxy population broadly covers the redshift range z ^ 0.03-2.6, and AGNs 
cover the redshift range of z % 0.1-8. These redshifts are based on both 
spectroscopic and photometric redshifts, and all sources at z > 4.76 have 
only photometric redshift estimates. We adopt the most probable photometric 
redshifts; however, some of the z > 4.76 photometric redshifts are consistent 
with being at lower redshifts (see Luo et al. 2010 for details). 

(A color version of this figure is available in the online journal.) 


models have been shown in Figures 5 and 6, respectively (con- 
tinuous curves). In Figure 7, we show sample probability dis- 
tributions for the best-fit dN/dS model parameters appropriate 
for AGNs detected in the SB. 

The bottom panels of Figures 5 and 6 show the fractional 
contribution that the respective cumulative and differential 
number counts of each source type make to the total counts. 
We find that AGNs still dominate the cumulative number 
counts at all flux levels in all four bands; however, at SB 
and FB fluxes below < 1(L 16 erg cm -2 s -1 , normal galaxies 
undergo a rapid increase in their numbers. At the faintest SB 
flux limit (»5.1 x 1 0 1 x erg cm~ 2 s -1 ), AGNs and normal 
galaxies reach source densities of ~ 1 4,900 and 12,700 deg -2 , 
respectively. These sky densities are factors of «^2.2 and 
~4.5 times larger than those reported at the ~2 Ms Chandra 
depth (see B04) for AGNs and normal galaxies, respectively. 
We find that our SB normal-galaxy number counts lie between 
the “optimistic” and “pessimistic” estimates from B04 and are 
in good agreement with previous work by Hornschemeier et al. 
(2003). Due to the small solid angle of our survey, the bright 
normal-galaxy number-count data drop off at SB fluxes >4 x 
1 0~ 1 6 erg cm 2 s _1 ; however, the extension of our best-fit power- 
law fit (red dotted curve in Figure 5(a)) agrees well with bright- 
end normal-galaxy number counts from the XMM-Newton 
serendipitous plus needles in the haystack surveys (Georgakakis 
et al. 2006) and the Chandra Multiwavelength Project (ChaMP; 
Kim et al. 2006). Normal galaxies constitute 46% ± 5% and 
43% ± 5% of the total cumulative number counts in the SB 
and FB, respectively, a result consistent with that found by 
XI 1, who showed that ?»40%-50% of the SB-detected sources 
within ~4' of the ~4 Ms CDF-S center are normal galaxies 
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Figure 5. Top panels of (a)-(d): cumulative number counts for the SB (a), HB (b), UHB (c), and FB (d) broken down into AGNs (open blue triangles), normal galaxies 
(open red squares), and Galactic stars (open green stars). The total number counts have been shown as filled circles. In each plot, the best-fit dN/dS parameterizations 
based on Equation (5) have been shown as blue long-dashed, red dotted, and green short-dashed curves for AGNs, normal galaxies, and Galactic stars, respectively. The 
total number-count model, based on summing the three contributing components, has been shown as a black curve. Bottom panels of (a)-(d): fractional contributions 
from AGNs, normal galaxies, and Galactic stars to the total number counts. For the majority of the flux ranges, AGNs dominate the number counts; however, normal 
galaxies provide significant contributions near the flux limits of the SB and FB. 

(A color version and a machine-readable table of the cumulative number-count data for the ~4 Ms CDF-S of this figure are available in the online journal.) 


(see their Figure 13(c)). This signifies more than a doubling in 
the fractional contribution that the normal-galaxy cumulative 
number counts make to the total number counts at the SB 
flux limit over that found in the Ms CDFs (see B04). In 
terms of the differential number counts, it appears that the 
contribution from normal galaxies exceeds that of AGNs at 
SB and FB fluxes below ^(1-2) x 1(1 17 erg cm 2 s -1 and 


^(3-7) xlO -17 erg cm -2 s _1 , respectively (see Figure 6(a)). 
We estimate that near the SB and FB flux limits, normal 
galaxies make up ~45%-70% of the differential number counts. 
The rapid increase in normal-galaxy counts signifies that deep 
Chandra surveys are on the verge of being normal-galaxy 
dominated at the faintest fluxes (see Section 5 for a discussion). 
In Table 2, we summarize properties of the number counts. 
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S = 4-8 keV Flux (ergs cm 2 s ') S = 0.5-8 keV Flux (ergs cm 2 s ') 

Figure 6. Differential number counts ( dN/dS ) vs. flux for the SB (a), HB (b), UHB (c), and FB (d), broken down into AGNs (open blue triangles), normal galaxies 
(open red squares), and Galactic stars (open green stars; for the SB. HB, and FB). The differential number counts have been estimated using flux bins of ~0.3 dex. In 
each plot, the best-fit dN/dS parameterizations based on Equation (5) have been shown as blue long-dashed, red dotted, and green short-dashed curves for AGNs, 
normal galaxies, and Galactic stars, respectively. The total number-count model, based on summing the three contributing components, has been shown as a black curve. 
Bottom panels of (a)-(d): fractional contributions from AGNs, normal galaxies, and Galactic stars to the differential number counts. In the SB and FB, we find that 
the differential number counts of normal galaxies appear to surpass those of AGNs for fluxes below ~(l-2) x !() 17 erg cm” 2 s~ 1 and ~(3-7) x 10 17 erg cm” 2 s , 
respectively. 

(A color version of this figure is available in the online journal.) 


In Figure 8, we show our derived total number-count 
distributions normalized by a Euclidean slope of S 1 5 
(i.e., N(> S ) x (SyS re f) 1 - 5 , where 5j-ef = 10 -14 erg cm -2 s -1 ). 
The Euclidean representation allows for comparisons with 
previous investigations at both the faint and bright ends of 


the distributions. For comparison, we have plotted number- 
count distributions from previous surveys with XMM-Newton 
(Brunner et al. 2008; Mateos et al. 2008; Cappelluti et al. 2009) 
and Chandra (Bauer et al. 2004; Kim et al. 2007; G08; Luo 
et al. 2008; Elvis et al. 2009; Puccetti et al. 2009). Some small 


8 


The Astrophysical Journal, 752:46 (23pp), 2012 June 10 


Lehmer et al. 



band number counts. These parameters have been summarized in Equation (5), and their values are provided in Table 1. 


Table 1 

Maximum-likelihood Best-lit Model Parameters 


Bandpass 


AGNs 



Normal Galaxies 

Galactic Stars 


IV- AGN 
A. i4 

oAGN 

Pi 

/i 2 agn 

/■AGN 

‘'break 

Kxg al 
A 14 


ex star 

A 14 

/3 star 

(i) 

(2) 

(3) 

(41 

(5) 

(6) 

(V) 

(8) 

(9) 

0.5-2 keV 

169.56 ± 8.69 

1.49 ±0.03 

2.48 ± 0.27 

6-0+Vi 

1.53 ± 0.10 

0 99 +O.O 8 
z - zz -0.07 

3.97 ± 0.37 

1 41 +O.W 
— 0.13 

2-8 keV 

573.13 ±27.49 

1.32 ±0.04 

9 cc+0.17 
-0.18 

6.4 ± 1.0 

1.10 ± 0.22 

2.29 ± 0.25 

0.64 ± 0.22 

1.79 ±0.50 

4-8 keV 

1463.27 ± 63.61 

1.08 ±0.10 

2.23 ± 0.09 

1 - 22 °o 7 8 

0.27 ±0.10 

2 40+0.38 


1 55 +019 
— 0.18 

0.5-8 keV 

562.20 ± 22.96 

1 04+0.04 

— 0.03 

2.35 ± 0.15 

81-i 5 4 

2.82 ± 0.26 

2 40+0. 1 1 
z>fU -0.12 

4.07 ±0.51 


Notes. Best-fit values and I rr errors for our number-count priors. Column 1 lists the bandpass. Columns 2-5 provide the double power-law differential 
number-count parameterization for AGNs, including the normalization A r ^ ( GN (Column 2; in units of 10 14 deg" 2 (erg cm~ 2 s -1 ) -1 ), faint-end slope /3 AGN 
(Column 3), bright-end slope /J AGN (Column 4), and break-flux / ( 4G ^ (Column 5; in units of 1 0 15 erg cm -2 s _1 ). Columns 6 and 7 provide the normal-galaxy 
single-power-law normalization A’jj' 1 and slope /3 gal , respectively. Columns 8 and 9 provide the Galactic star single-power-law normalization A - *™ and slope 
/3 star , respectively. These values were computed using the maximum-likelihood methods described in Section 2.2. 


differences in flux measurements (at the <10% level) are ex- 
pected due to both changes in calibration over the years for 
Chandra 15 and cross-calibration uncertainties between Chan- 
dra and XMM-Newton (see, e.g., Nevalainen et al. 2010). In 
general, we find good agreement between our number counts 
and those of previous studies with some minor differences. At 
faint SB fluxes, we find that our number-count measurements 
are in good agreement with those of G08 and B04, which are 
both based on the ~2 Ms CDF-N and ~ I Ms CDF-S; however, 
our estimates are »G0%-40% lower than those near the limit 
of the ~2 Ms CDF-S (L08). It is likely that this discrepancy is 
due to an incompleteness overcorrection of the faint-end counts 
in L08. In the HB, we find that our number counts are in good 
agreement with those of B04 but are somewhat lower than those 
of G08 and higher than those of L08 (at the ~ 1 0%-20% level 
in both cases). In the UHB, we find clear evidence of a break 
in the number counts at / 4 - 8 keV ~ 1 0 12 erg cnT 2 s _ 1 (see 
Figures 5 and 6). Previous investigations of the UHB in 
the CDF-S revealed no obvious break down to x 
10 -16 erg cm -2 s -1 (Rosati et al. 2002), which is very close to 
where we observe the break. We note that the G08 investigation 
found a break at a much brighter flux of ~6x 10 15 erg cm 2 s -1 . 

15 For example, as noted in XI 1, in 2009 January the Chandra ACIS-I 
ancillary response file (ARF) was updated, resulting in a flat ^9% reduction in 
effective area and a ~0%-8% reduction between 2 and 5 keV (see 
http://cxc.harvard.edu/ciao/why/caldb4. 1 . l_hrma.html). Further uncertainties, 
at the ~10% level, are expected due to spatial variations of contamination on 
the optical-blocking filter (see 

http://web.mit.edu/iachec/meetings/201 1/Presentations/Marshall.pdf). 


In the CDF-S, there are only «17 UHB sources brighter than 
this flux, which would make it difficult to identify such a break 
if present. In contrast, the G08 analysis made use of wide-area 
surveys like the EGS, ELAIS-N1, and XBOOTES and there- 
fore have better source statistics at bright UHB fluxes. In all 
four bands, the ~4 Ms CDF-S number-count measurements 
have large uncertainties at bright fluxes due to the relatively 
small solid angle of the survey; however, we find good general 
agreement between our measurements and those of previous 
investigations. 

In Sections 3.2 and 3.3 below, we further divide our number- 
count measurements into contributions from AGN and normal- 
galaxy subpopulations. For ease of presentation, we discuss only 
results from the SB and HB; however, full results from all four 
bandpasses are provided in the tables for reference. 

3.2. AGN Number Counts by Redshift, Intrinsic Column 
Density, and X-Ray Luminosity 

As discussed in Section 3.1, AGNs are the majority of the 
X-ray-detected sources and dominate the cumulative number 
counts over all fluxes. The ultra-deep Chandra data and redshift 
information allow for first-order estimates of the intrinsic 
absorption column densities, Ah, and 0.5-8 keV luminosities. 
Lx, of the detected AGNs. Values of Ah and Lx were estimated 
following the procedure outlined in Section 3.4 of Xue et al. 
(2010). The observed X-ray spectrum was modeled using xspec 
(Arnaud 1996) and an assumed absorbed power-law model 
(zpow x wabs x zwabs). The Galactic column density was 
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respectively, plus bright-end corrections from the Kim et al. (2007) number-count relations (see Section 3.4 for details). 
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S = 4-8 keV Flux (ergs cm 2 s *) S = 0.5-8 keV Flux (ergs cm 2 s *) 

Figure 8. Euclidean-normalized, cumulative number counts (N(> S) x (5/5 re f) 1-5 , where iS re f = 10 -14 erg cm -2 s _1 ) and 1 <t errors for the ~4 Ms CDF-S (black 
curves with gray error bars) for the SB (a), HB (b), UHB (c), and FB (d). The Miyaji & Griffiths (2002) predicted number-count boundaries (from ^1 Ms fluctuation 
analyses) are shown as the dotted regions in panels (a) and (b). Previous ^2 Ms CDF-N plus Ms CDF-S number-count estimates from B04 and ~2 Ms CDF-S 
number counts from L08 are indicated as blue dotted and red short-dashed curves, respectively. The G08 number counts, which include the ~2 Ms CDF-N and Ms 
CDF-S surveys, are indicated as dashed triple-dotted curves. Similarly, the ChaMP survey number counts from Kim et al. (2007), which include the ~2 Ms CDF-N 
and ~1 Ms CDF-S surveys, have been shown as green long-dashed curves. Number-count measurements from other XMM-Newton and Chandra surveys have been 
indicated, including the XMM-Newton serendipitous (Mateos et al. 2008; magenta triangles), Lockman hole (Brunner et al. 2008; brown squares), and COSMOS 
(Cappelluti et al. 2009; orange circles) surveys, as well as the Chandra COSMOS (Elvis et al. 2009; see also Puccetti et al. 2009; gold stars). In the case of the UHB, 
we highlight the ^1 Ms CDF-S number counts from Rosati et al. (2002; cyan dot-dashed curve). We note that all measurements shown (i.e., both curves and data 
points) correspond to measurements, and models are not plotted here. We find good overall agreement between our number-count measurements and those found in 
previous studies. 


set to 8.8 x 10 19 cm -2 (see Section 1), the intrinsic photon 
index was fixed at r; nt = 1.8 (e.g., Tozzi et al. 2006), and 
the best redshift estimate provided in the XI 1 catalog was 
adopted. Using this model, we determined the Ah value for 
each AGN that reproduced the observed ratio of count rates 
between the HB and SB. For sources detected only in the FB 
an observed spectral index of r o b s = 1.4 was adopted. This 
procedure will be effective for determining Ah and Lx values for 
AGNs in the Compton-thin/relatively unobscured regime (i.e., 
with Ah < 10 23 — 10 24 cm' 2 ) and is less reliable for a minority of 
sources that are expected to be Compton thick. In the Compton- 
thick regime, AGN spectra are significantly complex due to 
the combination of absorption, reflection, scattering, and line 


emission features that depend on metallicity, ionization state, 
and geometry (see, e.g., Murphy & Yaqoob 2009). Since most 
of our sources are in the low-count regime (i.e., <100 counts), 
such meaningful modeling is beyond our capabilities. Using 
our power-law model, estimated Ah, and intrinsic flux /x.int 
estimates, we computed Lx = 47r<i 2 /x i i n t(l + z) r ‘ nl_2 for 
each AGN that had a measured redshift (spectroscopic or 
photometric). 

In Figure 9, we display the breakdown of the AGN number 
counts in bins of redshift for the SB (left) and HB (right), 
with the middle panels showing the fractional contribution 
that sources in each redshift range make to the total AGN 
number counts. At the brightest fluxes, low-redshift {z < 1.5) 


11 


The Astrophysical Journal, 752:46 (23pp), 2012 June 10 


Lehmer et al. 



S = 0.5—2 keV Flux (ergs cm 2 s *) S = 2—8 keV (ergs cm 2 s x ) 

Figure 9. Top panels: cumulative AGN number counts (filled black circles) in the SB (left) and HB (right) broken down by contributions from AGNs in redshift 
ranges of z < 1.5 (open red circles), z = 1.5-3 (open green triangles), z > 3 (open blue squares), and no measured redshift (filled gray circles). AGN number-count 
predictions from the AGN population-synthesis models by G07 have been provided both with (dotted curves) and without (solid curves) an exponential decline in the 
luminosity function at z > 2.7. G07 AGN number-count models for all redshifts, z < 1.5, z = 1.5-3, and z > 3 have been shown as black, red, green, and blue curves, 
respectively. These models are discussed in detail in Section 4.1. Middle panels: fractional contribution that AGNs in each redshift range provide to the overall AGN 
number counts. The solid curves are the expectations from the G07 model without a declining XLF at z > 2.7. Bottom panels: the ratio of the AGN number-count 
data presented in the top panel divided by the G07 model (without a declining XLF at z > 2.7). 

(A color version and a machine-readable table of the redshift-divided AGN cumulative number-count data for all four bandpasses of this figure are available in the 
online journal.) 


AGNs dominate; however, going to fainter fluxes, we observe 
increasing contributions from z = 1.5-3 AGNs. At the flux 
limit of our survey, AGNs at z = 1.5-3 make up comparable 
contributions to the number counts as the z < 1.5 AGNs 
(i.e., ~40%). For AGNs at z > 3, we find a roughly steady 
fractional contribution to the number counts (^5%-15%) across 
the majority of the SB and HB flux ranges where they are 
detected. In Figure 10, we highlight the SB z > 3 AGN 
number counts and compare them with recent studies from 
the literature. We find that our number-count measurements 
for z > 3 AGNs are typically a factor of ~ 1 .5-2 higher than 
those found in the COSMOS survey fields (i.e., Brusa et al. 
2009; Civano et al. 2011) and in good agreement with those 
from the early Ms CDF-S study by Fiore et al. (2012). The 
difference between our values and those of the COSMOS studies 
is likely due to differences in z > 3 AGN selection techniques, 
variations in multiwavelength depth and photometric redshift 
completenesses, and field-to-field variance. 

In Figures 11 and 12, we show the breakdown of AGN 
number counts in bins of intrinsic absorption column density, 
/V[-| , and 0.5-8 keV luminosity. Lx, for the SB and HB. In both 
bandpasses, we find that unobscured X-ray-luminous sources 
dominate the bright-end number counts. Going to fainter fluxes, 
we find that more obscured and less X-ray-luminous sources 
make larger contributions. At the SB flux limit, AGNs with 
A h > 10"“ cm 2 and Lx < 1 0 43 erg s 1 make up the majority 
of contributions to the number counts. At the HB flux limit, 


we similarly find that AGNs with Ah > 10 22 cirT 2 dominate 
the number counts; however, in terms of luminosity, AGNs 
with Lx > 10 43 erg s~ and < 1 0 43 erg s -1 make comparable 
contributions, with the latter population increasing faster. 

3.3. Normal-galaxy Counts by Redshift and Morphology 

As noted in Section 3.1, normal galaxies compose a signifi- 
cant fraction («s32%^16%) of the source counts at the SB and 
FB flux limits. X-ray emission from normal galaxies is expected 
to be produced by a variety of populations that differ for dif- 
ferent galaxy populations (see, e.g., Fabbiano 1989, 2006 for a 
review). For late-type star-forming galaxies, populations asso- 
ciated with a combination of young and old stellar populations 
are expected to provide dominant contributions (i.e., high-mass 
and low-mass X-ray binaries, supernovae and their remnants, 
hot gas from starburst flows, and young stars; Colbert et al. 
2004; Iwasawa et al. 2009, 2011; Lehmer et al. 2010; Pereira- 
Santaella et al. 2011; Mineo et al. 2012). For passive early-type 
galaxies, X-ray emission is dominated by low-mass X-ray bina- 
ries (LMXBs) and hot X-ray-emitting gas (e.g., O’Sullivan et al. 
2001; Gilfanov 2004; Boroson et al. 201 1). In this section, we 
study the normal-galaxy number counts as a function of redshift 
and morphology. 

We measured the number counts of normal galaxies in two 
redshift bins divided at the median redshift Zmedian ~ 0.6. In 
Figure 13, we show the HB and SB number counts for 
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S — 0.5-2 keV (ergs cm 2 s *) 

Figure 10. Cumulative SB AGN number counts for sources at redshifts z > 3. 
The ^4 Ms CDF-S estimates and la error bars have been shown as open 
blue squares with error bars. If we add to the z > 3 AGN number counts the 
number counts for X-ray-detected sources without known redshifts, we obtain 
the open gray diamonds with error bars. Due to their lack of multiwavelength 
counterparts, in regions where extensive multiwavelength coverage is available, 
the sources with unknown redshifts are good candidates for being at z > 3. 
Recent constraints from Fiore et al. (2012) for the ~4 Ms CDF-S (filled black 
circles), Civano et al. (201 1) for Chandra COSMOS (filled green triangles), and 
Brusa et al. (2009) for XMM-Newton COSMOS (filled red squares) have been 
shown for comparison. The G07 models for z > 3 number counts have been 
shown for the cases where a high-redshift (z > 2.7) space-density decline is 
included (blue dotted curve) and is not (blue solid curve) included. The recent 
number-count predictions, based on the Aird et al. (2010) XLF measurement in 
the 2-10 keV band, are shown as a brown long-dashed curve. 

(A color version of this figure is available in the online journal.) 

normal galaxies at z < 0.6 and z > 0.6. In the SB, 
we find that low-redshift ( z < 0.6) normal galaxies dom- 
inate the bright-end counts; however, below flux levels of 
i=»(l-2) x UK 17 erg cm -2 s” 1 , high-redshift (z > 0.6) normal 
galaxies become the dominant galaxy population. By contrast, 
normal-galaxy number counts in the HB are dominated by low- 
redshift normal galaxies across the full flux range, with <20% 
of the number counts being from the high-redshift population. 

To characterize the number counts of normal galaxies divided 
by morphology, we made use of the multiwavelength data 
to classify the normal galaxies broadly as either late-type or 
early-type galaxies. We matched the 169 X-ray-detected normal 
galaxies to the optical source catalogs constructed by Xue et al. 
(2010), which provide rest-frame optical magnitudes based on 
SED fitting. To distinguish between relatively blue late-type 
galaxies and red early-type galaxies, we made use of the rest- 
frame U—V color. As described by Bell et al. (2004), the 
U and V bandpass pair straddle the 4000 A break and as a 
color provide a first-order indicator of the stellar age of the 
galactic stellar population. We applied the empirically calibrated 
redshift-dependent color division calculated in Section 5 of Bell 
et al. (2004) to separate active blue and passive red galaxy 
populations. This division was set at 

(U - V) rest = 1.15 - 0.3 lz - 0M(M v + 20.7), (8) 

such that galaxies redward and blueward of this division are can- 
didate early-type and late-type galaxies, respectively. We note 
that a single color division will not perfectly isolate the early- 


type and late-type galaxy populations. For example, edge-on or 
dusty late-type galaxies may be artificially reddened and there- 
fore classified as being early-type galaxies (see, e.g., Cardamone 
et al. 2010; Lusso et al. 2011). Furthermore, intrinsic scatter in 
the optical color distributions near the color cut will sometimes 
displace galaxies to the wrong side of the color division (e.g., 
“green” early-type galaxies that undergo low levels of recent 
star formation). We therefore visually inspected HST zsso-band 
images of all 169 galaxies to identify obvious cases where rest- 
frame colors provided incorrect morphological classifications. 
We found that 25 galaxies (»%5%) had obviously incorrect 
initial morphological classifications that were corrected. 

Out of the 169 normal galaxies, we classified 123 as late- 
type and 43 as early-type galaxies. In Figure 14, we show 
example HST three-color (B435, Woe, and z^so) images of 
10 early-type (top) and 10 late-type (bottom) galaxies that 
were in the normal-galaxy sample. Following the prescriptions 
discussed in Section 2, we computed number counts for each 
of these populations. In Figure 15, we show the number counts 
of the normal galaxies broken down by optical morphology. 
The normal-galaxy number counts for early-type and late- 
type galaxies are comparable at the brightest SB and HB flux 
levels. Progressing to lower fluxes, however, we find that the 
late-type galaxy population quickly rises and dominates the 
normal-galaxy number counts, and at the survey flux limits, 
late-type galaxies compose ~80% and »90% of the normal- 
galaxy number counts for the SB and HB, respectively (see 
lower panels of Figure 15). 

3.4. Contributions to the Cosmic X-Ray Background 

With the number-count estimates derived above, we 
can measure the corresponding contributions each source 
type and subtype make to the extragalactic CXRB. Here- 
after, we adopt CXRB intensities Dcxrb of (8.15 ± 
0.58) x 10” 12 erg cm' 2 s” 1 deg" 2 and (1.73 ± 0.23) x 
1 0“ 1 1 erg cm -2 s~ 1 deg” 2 for the SB and HB , respectively. These 
CXRB intensities were computed by summing components 
from (1) the ~ ] Ms CDF-S unresolved background intensity 
(D U n res; Hickox & Markevitch 2006), (2) the intensity of faint 
«1 Ms CDF-S sources (Qf a i n t; below s»5 x 10” 15 erg cm -2 s” 1 
and ~ 1 .4 x 10“ 14 erg cm 2 s” 1 for the SB and HB, respectively) 
as derived by Hickox & Markevitch (2006; see their Table 5), and 
(3) the bright-source intensity (Qbnght) derived from the best-fit 
number-count relations by Kim et al. (2007), which are based on 
~5500 X-ray sources from the ChaMP (e.g., Green et al. 2004; 
Kim et al. 2004a, 2004b). In the bright-source case, Dwight = 
yj>bnghi s'(dN/dS')dS', where .S'f a i nt is the brightest flux used to 

compute the faint-end counts, Slight = 10“ 11 erg cm -2 s” 1 (for 
both the SB and HB), and clN /dS’ is the best-fit differential 
number-count model from Table 3 of Kim et al. (2007). Thus, 
we derived Qcxrb = D unres + Dfaint + D hng ht for the SB and HB. 
For the FB and UHB, we adopted values of (2.54 ± 0.24) x 
10” 11 erg cm” 2 s” 1 deg” 2 (the sum of the SB and HB CXRB 
intensities) and (1.04 ± 0.23) x 10” 11 erg cm” 2 s” 1 deg” 2 , 
respectively. The UHB CXRB intensity was computed by con- 
verting the HB intensity to UHB assuming a power-law SED 
with T = 1.4 (see Moretti et al. 2009). 

In Table 2 (Columns 1 1—14), we provide the intensities 
and fractional CXRB contributions that sources in various 
subcategories in the ~4 Ms CDF-S provide. We find that the vast 
majority of the CXRB in each band is expected to be produced 
by AGNs, with normal galaxies producing only a small fraction 
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S = 0.5-2 keV Flux (ergs cm 2 s ') S = 2-8 keV Flux (ergs cm 2 s ') 

Figure 11. Same as Figure 9, but for AGNs divided into ranges of intrinsic column density (Jog( A'n/cm ~ 2 ) < 22 (open red circles), logtAn/cm 2 ) = 22-23 (open 
green triangles), log(/Y[[/cm 2 ) > 23 (open blue squares), and no measured Ah due to lack of redshift measurements (filled gray circles)). As in Figure 9, AGN 
number-count predictions from G07 have been provided as dotted and solid curves, corresponding respectively to the inclusion and exclusion of an exponential decline 
in the AGN space density at z > 2.7 (model curve colors correspond to symbol colors). 

(A color version and a machine-readable table of the An-divided AGN cumulative number-count data for all four bandpasses of this figure are available in the online 
journal.) 


(R^O.5%-3.0%) of the intensities. These results are similar to 
those found by B04, although revised to fainter flux levels. 

Due to the relatively small solid-angle coverage of the 
CDF-S survey («465 arcmin 2 ), we are unable to characterize 
the number counts for the relatively rare X-ray bright-source 
population (i.e., S > (1-10) x 10” 14 erg cm~ 2 s -1 ). To 
characterize the bright-end counts for each source type and 
subtype listed in Table 2 (Column 2) would require a wider 
survey with equivalent multiwavelength data to that used here 
for the ~4 Ms CDF-S. Despite this limitation, we can utilize 
our number-count estimates and the total bright-end CXRB 
intensities (i.e., Qbright) to estimate the total resolved CXRB 
intensities. In the first row of Table 2, we provide the bright-end- 
corrected resolved source intensities and fractional contributions 
to the CXRB. We find that «76%-88% of the CXRB intensity 
can be attributed to X-ray point sources, with the resolved 
fraction appearing to increase with median bandpass energy; 
however, a constant or decreasing resolved fraction with energy 
is not formally ruled out. Comparisons with past studies indicate 
a variety of levels of consistency with this result. For example, 
the studies of Hickox & Markevitch (2006) and Kim et al. 
(2007), which were used to estimate the total CXRB intensities 
and bright-end corrections for the CDF-S, find similar resolved 
CXRB fractions (using ^ 1-2 Ms depth CDF data) over the same 
energy ranges studied here. However, the studies of Moretti 
et al. (2003), B04, Worsley et al. (2005), and G08 find that 
the resolved CXRB fraction appears to decline somewhat with 
energy, using data that reach similar Chandra depths (» 1-2 Ms). 
These studies adopt results from either Moretti et al. (2003) 


and/or De Luca & Molendi (2004) when estimating the CXRB 
intensity and bright-end counts. We note, however, that even 
in these studies, the apparent decline in resolved CXRB with 
energy is less significant in the CDF-S itself (see, e.g.. Figure 2 of 
Worsley et al. 2005). We estimate that between the ~ 1-2 Ms and 
~4 Ms depths, the fraction of the resolved CXRB has increased 
by only «H%-2%, which is much smaller than the error bars 
of our measurements. Therefore, differences between previous 
resolved CXRB fractions are primarily related to either differing 
assumptions about the CXRB intensity and bright-end counts 
or field-to-field variations. In support of the latter point, L08 
find that the CDF-S number counts are lower than those of the 
CDF-N by ~25% : at the faintest fluxes. 

4. COMPARISONS WITH MODELS 

4.1. The AGN Number Counts and the Evolution 
of Accreting SMBHs 

The redshift-dependent AGN number counts (Figures 9 
and 10) provide a directly observable signature of the evolution 
of AGN activity in the universe. Since AGNs are the brightest 
X-ray-detected population, many previous investigations have 
focused on measuring their cosmic evolution (see Brandt & 
Hasinger 2005 for a review) and, in particular, the evolution of 
their X-ray luminosity function (XLF; e.g., Cowie et al. 2003; 
Ueda et al. 2003; Hasinger et al. 2005; Silverman et al. 2008; 
Ebrero et al. 2009; Aird et al. 2010). For example, Hasinger 
et al. (2005) made use of a variety of ROSAT, XMM-Newton, 
and Chandra surveys with highly complete optical/near- IR 
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S = 0.5-2 keV Flux (ergs cm 2 s ') S = 2-8 keV Flux (ergs cm 2 s ') 

Figure 12. Same as Figure 9, but for AGNs divided into ranges of intrinsic 0.5-8 keV luminosity (Jog( Lx/erg s ' 1 ) < 42 (open red circles), log(Lx/erg s _1 ) = 
42-43 (open green triangles), log(Lx/erg s” 1 ) = 43—44 (open blue squares), log(Lx/erg s ' 1 ) > 44 (open orange stars), and no measured Lx due to lack of redshift 
measurements (filled gray circles)). As in Figure 9, AGN number-count predictions from G07 have been provided as dotted and solid curves, corresponding respectively 
to the inclusion and exclusion of an exponential decline in the AGN space density at z > 2.7 (model curve colors correspond to symbol colors). 

(A color version and a machine-readable table of the Lx-divided AGN cumulative number-count data for all four bandpasses of this figure are available in the online 
journal.) 


redshift measurements to compute the observed evolution of the 
AGN XLF. This investigation found that X-ray-selected AGNs 
appear to undergo luminosity-dependent density evolution in 
their populations, where the space density of low-luminosity 
(Seyfert-type) AGNs peaks at z < 1 and the highly luminous 
(QSO-type) AGN space density peaks at z ^ 2. These direct 
measurements of the observed XLF evolution of AGN activity 
in the universe provide a useful constraint on the SMBFI 
accretion history; however, the observed X-ray bandpass used to 
make these measurements (i.e., ~0.5-2 keV) is susceptible to 
extinction, and therefore populations of AGNs with intrinsic 
X-ray absorption columns of Ah > 10 22 cm -2 will have 
underrepresented intrinsic accretion activity. 

Using the unresolved hard X-ray (~3— 1 00 keV) background 
intensity spectrum as a constraint on the obscured (Ah ^ 
10 22 -10 24 cm" 2 ) and Compton-thick (A H > 10 24 cm" 2 ) AGN 
populations, Gilli et al. (2007; hereafter G07) used the observed 
XLF evolution from Hasinger et al. (2005), local distributions 
of the luminosity-dependent obscured-to-unobscured AGN 
fraction, and a model distribution of AGN population column 
densities to form phenomenological models describing the in- 
trinsic evolution of the AGN XLF (see also Ueda et al. 2003; 
Draper & Ballantyne 2009; Treister et al. 2009 for additional 
prescriptions). These models make direct predictions for the 
number-count distributions of sources with different ranges 
of z, Ah, and intrinsic 0.5-8 keV luminosity Lx 16 and can 


16 The G07 phenomenological model number-count predictions can be found 
at http://www.bo.astro.it/~gilli/counts.html. 


be directly compared with our observed AGN number-count 
subpopulations (see Figures 9-12). 

In Figures 9-12, we have highlighted the number-count 
predictions for z. Ah, and Lx selected AGNs by G07 appropriate 
for the Hasinger et al. (2005) XLF as observed (solid curves) 
and after applying a high-redshift ( z > 2.7) exponential 
decline to the XLF (dotted curves; see Schmidt et al. 1995 
for motivation). In general, we found better agreement by not 
including the exponential declining term (i.e., the solid curves in 
Figures 9-12). However, since the majority of the AGNs in this 
study are at z < 2.7, there is generally little difference between 
the models that include and do not include an exponentially 
declining term. For clarity in comparisons, in Figures 9, 11, 
and 1 2 we have provided lower panels showing the ratio between 
our observed number counts and the G07 models without the 
exponential declining term. 

For redshift-selected number counts, we found reasonable 
agreement (within a factor of ~2) between the observed number 
counts and the G07 model across the entire SB and HB flux 
ranges (see Figure 9). We note that a non-negligible fraction of 
the AGN population does not have redshifts available. Since the 
multiwavelength coverage over the entire CDF-S is extensive 
and deep, sources without redshifts (neither spectroscopic nor 
photometric) have very faint optical/near-IR counterparts and 
are therefore good candidates for high-redshift (z > 3) AGNs. 
If we add the number counts from AGNs without redshifts to the 
Z > 3 number counts, we find improved agreement with the G07 
model prediction (see Figure 10, diamonds). For comparison, we 
have also plotted the predictions from the recent work by Aird 
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S = 0.5-2 keV Flux (ergs cm 2 s ') S = 2-8 keV Flux (ergs cm 2 s ’) 

Figure 13. Top panels: cumulative normal-galaxy number counts (filled black circles) in the SB (left) and HB (right) broken down by contributions from galaxies at 
redshifts z < 0.6 (open orange stars) and z > 0.6 (open green upside-down triangles). The redshift division between the two normal-galaxy samples was made at the 
median galaxy redshift. The dotted curves in the left panel show the predicted z < 0.6 (orange dotted curve) and z > 0.6 (green dotted curve) number counts based 
on the observed XLF pure-luminosity evolution parameterization derived by Ptak et al. (2007) using observed XLFs out to z ~ 0.75. At i > 0.6, these predictions 
are based on an extrapolation of the XLF parameterization out to z ** 2. The solid curves in the left and right panels show the predicted number counts for - < 0.6 
(orange solid curves) and z > 0.6 (green solid curves) galaxy populations based on the observed evolution of stellar-mass functions converted to XLFs using X-ray 
scaling relations (see Section 4.2 for details). Middle panels: the fractional contributions that z < 0.6 and z > 0.6 galaxy populations make to the total galaxy number 
counts. Bottom panels: the ratio between the data and predictions based on the solid curves in the top panels. 

(A color version and a machine-readable table of the redshift-divided normal-galaxy cumulative number-count data for all four bandpasses of this figure are available 
in the online journal.) 


et al. (2010), which is based on 2-10 keV XLF measurements 
from the CDFs and the AEGIS-X surveys. The Aird et al. (2010) 
predictions are also consistent with our z > 3 AGN number 
counts; however, if all of the AGNs with unknown redshifts are 
in the z > 3 range, our number counts become more consistent 
with the G07 model. Regardless of the redshifts of the sources 
lacking redshift identifications, our number counts appear to be 
a factor of >2 higher than those predicted by the G07 model 
with exponential decline. 

For the number-count predictions based on Ah, we find 
significant discrepancies (factors of ss 1.5-3) between those 
observed and those predicted by the G07 model (see bottom 
panels of Figure 11). This situation is similar regardless of 
whether the exponential decline is applied or not. We suspect 
that some of the disagreement here will be due to (1) cosmic 
variance of the populations, (2) our simplistic calculation of 
Ah based on band ratios (see discussion in Section 3.2), which 
differs from the spectral models from G07 (see their Figure 1), 
(3) the distribution of Ah values for sources without redshifts, 
and (4) imperfections in the G07 model descriptions. Regarding 
point (3), we note that sources with unknown redshifts have 
a median effective photon index of r me dian ~ 0.7, which is 
somewhat harder (softer) than the r me( i; an 1 .0 (T me( ji an 0.4) 
for sources with log(A H /cm _2 ) = 22-23 (log(A H /cm _2 ) > 23). 
Therefore, we expect that most of the sources with unknown 
redshifts will have log(An/cm " 2 ) > 23 regardless of their 


redshifts. This would improve the agreement between observed 
and model number counts for log(/V H /cm -2 ) > 23. Despite 
the apparent disagreements and limitations, we do find basic 
agreement in the general trend: more heavily obscured AGNs 
have larger contributions at fainter fluxes. 

For number counts divided by intrinsic 0.5-8 keV luminosity 
Lx, we find good (within a factor of ~2) basic agreement 
between our measurements and the G07 model (see Figure 12); 
however, this is not true for the decline model (see dotted curves 
in upper panels of Figure 12). 

In summary, we find that the G07 phenomenological models 
provide a good description of the redshift, column density, 
and intrinsic luminosity distributions to the extent that we can 
constrain these values observationally. Better measurements of 
the column densities and Lx would be needed to improve these 
models further. However, unless a significant fraction of objects 
with photometric redshifts z > 3 prove to be lower-redshift 
contaminants, we do find that applying an exponential decline 
to the AGN XLFs at z > 2.7, under the G07 parameterizations, 
provides a poor characterization of the redshift- and Lx-divided 
number counts. 

4.2. The Rapidly Increasing Normal-galaxy Counts 

In Section 3.3, we presented number-count estimates for 
the normal-galaxy population and subpopulations separated by 
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Figure 14. Three-color (£ 435 , Vm, and 3350 ) HST images of a sample of X-ray-detected normal galaxies ordered by redshift. The top panels show 10 example late-type 
galaxies, and the bottom panels show 10 example early-type galaxies. The angular extent of each image has been chosen to subtend a physical distance of «j 30 kpc on 
a side in the rest frame of each galaxy. The redshift of each source is indicated in the upper left, and a 1 arcsec bar has been placed in the lower-left hand comer to show 
angular scale. The images are centered on the coordinates corresponding to each galaxy, which were taken from Columns 1 8 and 1 9 of Table 3 in X 1 1 . These positions 
were obtained primarily from VLT WFI observations, and therefore small (<0. 1-0.2 arcsec) offsets are sometimes seen in these HST images (note the angular scale). 


redshift and galaxy morphology. Given that these X-ray source 
populations are uniquely accessible to the ultra-deep survey 
regime probed by the Chandra Deep Fields, the evolution of 
the normal-galaxy XLFs has not been as well characterized 
as those of AGN populations (see Section 4.1). Flowever, 
initial investigations have placed first-order constraints on the 
evolution of normal-galaxy XLFs out to z ^ 0.75 (e.g., Norman 
et al. 2004; Georgakakis et al. 2007; Ptak et al. 2007; Tzanavaris 
& Georgantopoulos 2008). These investigations have provided 
evidence for a rapidly increasing late-type galaxy 0.5-2 keV 
XLF that is consistent with pure-luminosity evolution, where 
L* oc (1 + z) h , and b 1.5-3. By contrast, the early-type 
galaxy 0.5-2 keV XLF appears to be changing more slowly 
with redshift, and in the Tzanavaris & Georgantopoulos (2008) 
study, the XLF is consistent with being constant or decreasing 
with redshift (the value of b ranges from ~ —2.3 to + 2.6 in 
these studies). However, these results are based on shallower 
data than the ~4 Ms CDF-S and thus far provide only estimates 
for how the bright end of the galaxy XLFs evolved out to 
z ~ 0.75. In the ~4 Ms CDF-S, ~29% of the normal galaxies 
are at z > 0.75, and the majority of these sources are only 
accessible by this survey. In the paragraphs below, we interpret 
the number counts of the normal-galaxy subpopulations (e.g., 
see Figures 13 and 15) in the context of what we expect 
from previous estimates of the galaxy XLF evolution and the 


evolution of galaxy properties (e.g., SFR and stellar mass) and 
X-ray scaling relations (e.g., the X-ray/SFR relation). 

Following the formalism of Ranalli et al. (2005), the normal- 
galaxy number counts can be expressed as the following double 
integral: 

r Zmax nL max V 

N(>S) = X deg / dz (p x ( log L x , z) d log L x ——, 
JZMn Jl^s) dzdQ. 

1 (9) 

where X s d r eg = 3.05 x 1 0 4 converts from sr” 1 to deg” 2 , < px 
is the redshift-dependent XLF, dV/dzdU. is the cosmology- 
dependent differential volume element (comoving volume per 
unit redshift per unit solid angle), and [z m j n , z max ] and [L m i n (5'), 
L max ] are, respectively, redshift and luminosity integration lim- 
its. The value of L m i n (S) is the rest-frame luminosity corre- 
sponding to a flux S and redshift z . In the upper-left panels of Fig- 
ures 13 and 15, the dotted curves show 0.5-2 keV number-count 
predictions based on previous measurements of the evolution of 
galaxy XLFs for redshift- and morphology-divided galaxy popu- 
lations, respectively. These curves were computed using the Ptak 
et al. (2007) redshift-dependent Schechter parameterizations of 
the observed galaxy XLFs, which assume pure-luminosity evo- 
lution of the form L* oc (1 + z ) 2 ' 3 and (1 + z ) L6 for late-types 
and early-types, respectively. We also experimented with using 
a log-normal parameterization of the galaxy XLF and found 
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S = 0.5-2 keV Flux (ergs cm 2 s ! ) S = 2-8 keV Flux (ergs cm 2 s *) 

Figure 15. Top panels: similar to Figure 13, but for normal galaxies divided into late-type (open blue squares) and early-type (open red triangles) galaxy populations. 
The classification of galaxy morphology is explained in Section 3.3. The dotted curves in the upper-left panel show the predicted late-type (blue dotted curve) and 
early-type (red dotted curve) galaxy number counts based on the observed XLF pure-luminosity evolution parameterization derived by Ptak et al. (2007) using observed 
XLFs out to z ^ 0.75. These predictions are based on extrapolation of the XLF parameterization out to z ^ 2. The solid curves in the left and right panels show the 
predicted number counts for late-type (blue solid curve) and early-type (red solid curve) galaxy populations based on the observed mass functions converted to XLFs 
using X-ray scaling relations (see Section 4.2 for details). Middle panels: the fractional contributions that late-type and early-type galaxies make to the total galaxy 
population. Bottom panels: the ratio between the data and predictions based on the solid curves in the top panels. The reasonable agreement between the number-count 
data and models, both here and in Figure 13, indicates that the rapidly rising galaxy counts are largely attributed to the evolution of galaxy physical properties (e.g., 
stellar mass and star formation rate) with little redshift evolution in scaling relations between these physical properties and X-ray power output. 

(A color version and a machine-readable table of the morphology-divided normal-galaxy cumulative number-count data for all four bandpasses of this figure are 
available in the online journal.) 


similar results. When estimating galaxy number-count predic- 
tions, we used Equation (9) and extrapolated the Ptak et al. 
(2007) parameterizations out to z ~ 2. The Ptak et al. (2007) 
XLFs were directly measured for galaxies out to z ~ 0.75 (using 
galaxies with redshifts as large as z ^ 1 .2) and are in good agree- 
ment with other galaxy XLF measurements from the literature 
(e.g., Norman et al. 2004; Tzanavaris & Georgantopoulos 2008). 
The uncertainties on the normalization of the XLF parameteriza- 
tion are a factor of ~2. We find that, within the uncertainties, the 
XLFs provide reasonable predictions for the z < 0.6 (Figure 13, 
orange dotted curve and stars) and late-type galaxy (Figure 15, 
blue dotted curve and squares) number counts. For these two 
samples, the largest differences come from the z < 0.6 sample 
at 0.5-2 keV fluxes >3 x 10' 17 erg cm -2 s -1 , which may be 
due to differences in galaxy classification of X-ray-detected 
sources, cosmic variance, and statistical errors associated with 
the XLF fits. We find that the extrapolated XLFs appear to 
overpredict the z > 0.6 (Figure 13, green dotted curve and 
upside-down triangles) and early-type galaxy (Figure 15, red 
dotted curves and triangles) number counts by more than a 
factor of two, suggesting that our extrapolation of the early- 
type galaxy XLF evolution out to z ^ 2 is an overestimate. 
Therefore, it is possible that the rate by which the early-type 
galaxy XLF rises with redshift decreases somewhat going from 


Z ss 0.75 to 2 and/or the adopted rate of increase is too high (e.g., 
see Tzanavaris & Georgantopoulos 2008, who find a lower rate 
of increase). To obtain independent predictions of the galaxy 
XLFs and number counts for galaxies out to z rs 2, we can 
use the combination of X-ray scaling relations (e.g., the 
X-ray/SFR relation) and the redshift evolution of galaxy phys- 
ical properties (e.g., SFR and stellar mass) that are measured 
using multiwavelength observations. 

Recent X-ray stacking investigations of large normal-galaxy 
populations selected by morphology and galaxy physical prop- 
erties (e.g., SFR and stellar mass) have now provided insight into 
how the X-ray emission from normal galaxies evolves relative to 
galaxy physical properties out to z 1-4. Lehmer et al. (2008) 
reported that late-type star-forming galaxies, the most numerous 
galaxy population in the universe, have mean Lx/SFR values 
that are roughly constant out to z « 1 .4. Further X-ray stacking 
of distant z ~ 1.5-4 Lyman break galaxies appears to show 
similar mean Lx/SFR values (e.g., Brandt et al. 2001; Lehmer 
et al. 2005; Laird et al. 2006; Cowie et al. 2012; A. Basu-Zych 
et al. 2012, in preparation). Similarly, Vattakunnel et al. (2012) 
have studied z < 1.2 normal galaxies in the ~4 Ms CDF-S 
at X-ray and radio wavelengths and found that the Lx-L 1.4 ghz 
correlation for these galaxies is consistent with that expected 
from the local Lx/SFR relation, albeit with significant scatter 
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(~0.4 dex). Additional indirect support for a roughly constant 
L x /SI R ratio with redshift has been provided by Dijkstra et al. 
(2012), who used the observed SFR density evolution of the 
universe to show that the unresolved SB CXRB can be fully 
explained by an Lx/SFR ratio oc(l + z) b , where b is constrained 
to be less than 1 .4. 

For early-type galaxies, Lehmer et al. (2007) and Danielson 
et al. (2012) showed that the SB X-ray luminosity per unit 
/1-band luminosity (Lx/ L#) appears to undergo only mild 
evolution (a(l + z ) L2 ) to z « 1.2 for the most optically 
luminous galaxies, which is expected to be due to X-ray- 
emitting hot gas. In contrast, however, lower-luminosity early- 
type galaxies have been found to have rising mean Lx/ Lb 
values with increasing redshift, which is expected to be due to a 
fading LMXB population toward the present day (Lehmer et al. 
2007; A. E. Hornschemeier et al. 2012, in preparation). Similar 
expectations have been found from X-ray binary population- 
synthesis models, which predict an evolution of Lx/ L# oc 
(1 + z) 2 at z < 1 (Fragos et al. 2012, in preparation; see also 
Ghosh & White 2001). 

The above investigations provide estimates for how scaling re- 
lations between X-ray luminosity and galaxy physical properties 
have evolved out to z ^ 1-4, where the majority of the ~4 Ms 
CDF-S normal galaxies are detected. It seems that only the 
scaling relation between LMXB emission and stellar mass un- 
dergoes significant evolution with cosmic time; however, these 
populations are expected to dominate only in early-type galax- 
ies, which make up a minority of the X-ray-detected sources 
in the ~4 Ms CDF-S. It is therefore likely that the evolving 
normal-galaxy XLF and rising galaxy number counts, observed 
here and in previous studies, are largely due to the rapid evo- 
lution of galaxy physical properties. Recent multiwavelength 
surveys (e.g., COSMOS; Scoville et al. 2007) have provided 
significant new constraints on the evolution of optical/near-IR 
luminosity, SFR, and stellar-mass functions. These measure- 
ments can be combined with scaling relations between X-ray 
luminosity and galaxy physical properties to make predictions 
for the evolution of the normal-galaxy XLFs and number counts 
for galaxies out to z ^ 2. 

Following Avni & Tananbaum (1986), (px can be estimated 
through the following transformation: 

/ OO 

<A>( log &, z)P (log L x | log &) d log iL, 

-OO 

( 10 ) 

where i) represents a galaxy physical property (e.g., optical 
luminosity, stellar mass, or SFR) by which ^(log i), z) has been 
measured, and LTlog L x | log d) is the probability distribution 
for observing Lx given the value t) . 

We make use of the redshift-dependent stellar-mass functions 
measured by Ilbert et al. (2010) for late-type and early-type 
galaxies out to z ^ 2 (i.e., f) — AL*). To convert < 7 >M„(log AL*, z) 
to ^x(logLx,z) via Equation (9) requires knowledge of the 
transformation between X-ray luminosity and stellar mass (i.e., 
L’(log Lx| log M^). This function is expected to depend on 
galaxy morphology. In its most general form, we can express 
the transformation function as follows: 


P (log L x I log M*) 


x/2tzo 


x exp 


(log L X ,mod(ALT - log Lx) 2 ' 


2cr 2 


( 11 ) 


where L x , mo d(AL*) is the predicted X-ray luminosity given a 
value of the stellar mass AL* and a represents the scatter in 
the relation. For normal late-type galaxies in the local universe, 
the X-ray luminosity has been shown to correlate strongly with 
SFR (e.g., Persic & Rephaeli 2007; Lehmer et al. 2010; Pereira- 
Santaella et al. 20 1 1 ; Symeonidis et al. 20 1 1 ; Mineo et al. 20 1 2), 
and as discussed above, this correlation is observed to hold out 
to z > 1. We made use of the Lehmer et al. (2010) X-ray /SFR 
relation for late-type star-forming galaxies in the local universe: 

L 2 — iokeV = aAL* + /3SFR, (12) 

where a = 9.08 x 10 28 erg s -1 ALq 1 and fi = 1.62 x 10 39 erg s -1 
(AL G yr -1 ) -1 , AL, is in units of AL G , SFR is in units of Af 0 yr -1 , 
and L 2 — iokeV is in units of erg s . Studies of local and distant 
galaxy populations have shown that the SFR for late-type star- 
forming galaxies is further strongly correlated with AL* out to 
Z ^ 2 (e.g., Daddi et al. 2007; Elbaz et al. 2007; Peng et al. 
2010), albeit with a rapidly evolving normalization (i.e., mean 
SFR/AL*) with redshift (e.g., Zheng et al. 2007; Karim et al. 
2011; Cen 2011). To convert SFR to AL, for late-type star- 
forming galaxies, we adopted the relation provided by Karim 
et al. (2011): 

SFR/AL* = c 0 (l +z) 3 ' 5 (AL*/10 n AL G / SFR , (13) 

where c 0 ~ 0.0263 Gyr -1 and /3sfr ~ —0.4 are fitting 
constants. Combining this relation with Equation (12), we arrive 
at the following expression: 


L 2 — iokeV = 0 / AL* + y( 1 + z) 35 AL*- 6 , (14) 

where y = 1.07 x 10 33 , AL* is in units of AL 0 , and L 2 -iokeV is 
in units of erg s _1 . When computing number-count predictions, 
we converted the 2-10 keV luminosities provided here to other 
bandpasses assuming a power-law SED with T = 1.9, the 
average stacked SED of the X-ray-detected normal galaxies 
(Young et al. 2012). The overall scatter (accounting for both 
SFR/AL* and L X /SFR relations) in this relation is estimated to 
be a zz 0.4 dex. 

For normal early-type galaxies, the X-ray luminosity has been 
shown to correlate with A'- band luminosity, L K , which provides 
a direct proxy for galaxy stellar mass (e.g., Gilfanov 2004; 
Boroson et al. 2011). We converted the Boroson et al. (2011) 
Lx/Lx relations for early-type galaxies to Lx/AL* relations 
assuming a single mass-to-light ratio logAL*/Lx ~ —0.2 
(where AL* and Lk are in solar units) characteristic of early- 
type galaxies. These relations are separated into contributions 
from LMXBs and hot X-ray-emitting gas: 

log L 0 .3-8kev(LMXB) 29.2 + log AL* + 2 log(l + z) 
logL 0 . 3 - 8ke v(gas) ~ -6.07 + 4.03 log AL*, (15) 

where the 21og(l +z) term accounts for the expected redshift 
evolution for LMXBs in early-type galaxies (see discussion 
above). When computing number-count predictions, we con- 
verted the 0.3-8 keV relations to other bandpasses assum- 
ing a power-law SED with T = 1.8 for LMXBs and an 
apec plasma model with an Lo. 3 - 8 keV-dependent temperature 
(L x oc Lj]^ skeV ) and solar abundances for the hot gas. The 
scatter in the LMXB and hot-gas relations is estimated to be 
R^O.2 dex and ~ 1 .0 dex, respectively (see Boroson et al. 2011). 

The relations expressed in Equations (14) and (15), combined 
with the stellar- mass functions from Ilbert et al. (2010), were 
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folded into Equation (9) to obtain redshift-dependent XLFs for 
late-type and early-type galaxies, respectively. We measured 
the number counts of both populations in two redshift bins by 
integrating Equation (9) over the redshift intervals z < 0.6 
and z > 0.6 and summing the contributions from each galaxy 
morphology type. In Figures 13 and 15, we show our model 
number counts in bins of redshift and morphology, respectively, 
and in the bottom panels, we plot the ratio of our data to the 
model. We remind the reader that the models implemented here 
are not fits to the data in any way. For the redshift-divided 
normal-galaxy samples, we find good agreement between our 
data and the models at the SB and HB flux limits; however, 
for the number counts for normal galaxies at z > 0.6, we 
find significant differences at the bright end (/o.5-2keV and 
/2-8keV > 10 -16 erg crrT 2 s" 1 ). The bright-end number counts 
for z > 0.6 galaxies are based on only ~ 1-5 objects; therefore, 
the differences with the models may be due to either cosmic 
variance or misclassification of some fraction of these sources 
(i.e., if some of these normal galaxies are actually AGNs). 
Furthermore, the general trends seen in the data and models (e.g., 
the increasing contributions and faint-end dominance in the SB 
of z > 0.6 galaxies) appear to be broadly consistent with each 
other. For the morphology-selected number counts, we find very 
good agreement (less than a factor of « 1 .5 difference) between 
the data and the model with the exception of the HB early- 
type galaxy counts. We note, however, that the HB number- 
count measurements for early-type galaxies are based on small 
numbers of galaxies (four galaxies) and are sensitive to the 
classifications discussed in Section 3.1 above. 

In summary, we find that our models provide good charac- 
terizations of the X-ray number counts in redshift bins for both 
late-type and early-type galaxy populations down to the flux 
limits of our survey (with the exception of z > 0.6 normal 
galaxies with bright SB and HB fluxes and HB early-type galaxy 
counts; however, see discussion above). For late-type galaxies, 
this result suggests that X-ray scaling relations in the local uni- 
verse appear to hold out to at least z ~ 1-2, in agreement with 
previous investigations. For early-type galaxies, these results are 
consistent with there being little redshift evolution in how hot- 
gas X-ray emission scales with stellar mass and modest redshift 
evolution (oc(l + z) 2 ) in the LMXB luminosity per stellar mass. 

5. THE IMPENDING DOMINANCE OF NORMAL 
GALAXIES: PREDICTIONS FOR FUTURE DEEPER 
X-RAY OBSERVATIONS 

The number-count results presented above show that normal 
galaxies play an increasingly dominant role in the extragalactic 
X-ray source population going to fainter fluxes. Our observa- 
tions show that the cumulative SB and FB number counts of nor- 
mal galaxies are on the brink of surpassing those of AGNs (see 
Figure 5), and the slope of the galaxy number counts ( dN/dS ) 
already exceeds that of AGNs (see Figure 6). At the SB flux 
limit of the ~4 Ms CDF-S normal galaxies are reaching the 
unprecedented source density of ~ 1 2,700 deg -2 , which consti- 
tutes ~46% of the total number counts. Given the comparable 
source densities of normal galaxies and AGNs at the survey flux 
limit and the relatively sharp faint-end rise in the galaxy number 
counts versus AGNs ( N(> S) power-law slope of /( gal — 1 ~ 1.2 
versus /i | ' GX — 1 ~ 0.5; see Table 1), it is clear that going to 
fainter Chandra fluxes beyond ~4 Ms will yield larger gains 
in X-ray source densities than previous increases in depth. XI 1 
showed that exposures of the CDF-S region will yield substantial 
gains in sensitivity. For example, a feasible Ms exposure 



S = 0.5-2 keV Flux (ergs cm 2 s *) 

Figure 16. Measured SB number counts (filled black circles) with contributions 
from AGNs (open blue triangles) and normal galaxies (open red squares) shown. 
Our best-fit parameterizations from Equation (5) for all sources, AGNs, and 
normal galaxies have been shown as solid black, dashed blue, and dotted 
red curves, respectively. An extrapolation of these models indicates that the 
normal-galaxy number counts will overtake those of AGNs at SB flux levels of 
f=a(3— 5) x KT 18 erg cm -2 s" 1 . The ultimate flux limits of the ~4 Ms CDF-S, 
a potential ~10 Ms CDF-S, and future ?»4 Ms SMART-X exposures have been 
indicated. These extrapolations indicate that a «10 Ms CDF-S would have SB 
number counts at the flux limit that reach source densities of ~46,000 deg~ 2 
and are dominated by normal galaxies. More ambitiously, a ^4 Ms SMART-X 
exposure would reach source densities of «;500,000 deg 2 and would have a 
^90% contribution from normal galaxies. For reference, we have plotted the 
source densities of galaxies with R-band magnitudes brighter than R ^ 24 and 
R Ri 27 (gray short-dashed lines). 

(A color version of this figure is available in the online journal.) 


would permit a factor of «H,8 increase in depth (in terms of 
flux limit) over the current ~4 Ms exposure. Therefore, explor- 
ing the normal-galaxy-dominated regime of the number-count 
distribution is a sensible goal that may be pursued in future 
surveys. At present, the only X-ray mission capable of reaching 
these fluxes is Chandra (see below), and this will likely be the 
case for at least the next two decades. 

In Figure 16, we provide an expanded view of the SB number 
counts and their contributions by extragalactic source type 
(AGNs and normal galaxies) in the faint-flux regime. Our best- 
fit number-count parameterizations (based on Equation (5)) 
have been indicated. Direct extrapolation of these parameter- 
izations indicates that the number counts of normal galaxies 
will overtake those of AGNs just below our current flux limit 
at «(3-5) x 10 -18 erg cm -2 s -1 . With a Ms exposure in 
the CDF-S, we could study number counts down to a sensi- 
tivity limit of «i2.8 x 10~ 18 erg cm -2 s -1 (see vertical dotted 
line in Figure 16), although we expect that the faintest sources 
will have most probable fluxes of «^5.2 x 10~ 18 erg cm -2 s -1 
(see Section 2.1 for distinction). At such a depth, we would ex- 
pect a factor of « 1.7 increase in total X-ray source densities to 
~46,000 deg -2 (based on extrapolation of our best-fit number- 
count parameterization), with normal galaxies making up «57% 
of the limiting source density. After considering variations in 
sensitivity across the field, we expect that a 10 Ms exposure 
would yield ~ 1 020- 1 080 total X-ray-detected sources. 
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While Chandra is the only operating observatory capable of 
observing to the ultra-deep regime studied here, the realiza- 
tion of future X-ray mission concepts like Generation-X and 
SMART-X , 1 7 which have much larger light-collecting areas and 
comparable or improved imaging resolution over Chandra, 
would allow for the exploration of new populations of extra- 
galactic X-ray sources, including new populations of distant 
(z > 0.6) normal galaxies (see, e.g., Figure 13). Under its 
current specification, a ~4 Ms SMART-X survey would detect 
SB sources with fluxes as low as ~3 x 10~ 19 erg cm' 2 s' 1 . 
Direct extrapolation of our number counts to these levels pre- 
dicts source densities of ^500,000 deg' 2 , with ~90% contribu- 
tion from normal-galaxy populations. 

6. SUMMARY 

In this paper, we have provided new quantitative measure- 
ments of the X-ray number counts in the ~4 Ms CDF-S using 
four bandpasses: 0.5-2 keV, 2-8 keV, 4-8 keV, and 0.5-8 keV 
(SB, HB, UHB, andFB, respectively ). Our analyses focus on the 
faintest flux regimes only accessible by such ultra-deep Chandra 
surveys. We draw the following key conclusions: 

1. We make use of a Bayesian approach, flux probability dis- 
tributions, and maximum-likelihood techniques to measure 
number counts for sources detected in the ~4 Ms CDF-S 
down to SB, HB, UHB, and FB flux limits of 5.1 x HU 18 , 
3.7 x 10' 17 , 4.6 x 10~ 17 , and 2.4 x 10' 17 erg cm' 2 s _1 , re- 
spectively; these are factors of « 1.9-4. 3 times fainter than 
previous number-count investigations. 

2. At the flux limits of our survey, we reach source densities 
of 27,800, 10,500, 8400, and 22,600 deg' 2 in the SB, 
HB, UHB, and FB, respectively. AGNs are the majority 
contributors to the total number counts at all fluxes in all 
four bands (especially in the HB and UHB). In the SB, we 
reach AGN sky densities of R= 1 4,900 deg' 2 , the highest 
reliable AGN sky density measured at any wavelength. 

3. The normal-galaxy number counts rapidly rise, compared 
with AGNs, and at the SB and FB flux limits, normal 
galaxies make up «46% and ~43% of the number counts, 
respectively. At the limiting SB flux, normal galaxies reach 
a sky density of ~ 1 2,700 deg' 2 , which is a factor of 
«4.5 times higher than measured in the ~2 Ms surveys. 

4. Stated in terms of the intensity on the sky. X-ray-detected 
sources in the ~4 Ms CDF-S can account for ~76%-88% 
of the CXRB. We estimate that the increased exposure from 
Rs 1-2 Ms to ~4 Ms has resulted in an increase of ~ 1 %—2% 
resolved CXRB. 

5. In the SB and HB, AGNs at high redshifts (z > 1.5), with 
moderate to significant obscuration ( /V H > 10 22 cm' 2 ), 
and relatively low intrinsic 0.5-8 keV luminosity (Lx < 
10 43 erg s' 1 ) make increasingly important contributions to 
the AGN number counts going to fainter fluxes, and at the 
flux limits, these sources make up >40% of the overall 
AGN number counts. These trends are broadly consistent 
with those expected from the phenomenological models 
of G07. 

6. In the SB and HB, late-type star-forming galaxies with 
Z < 0.6 provide majority contributions to the normal- 
galaxy number counts, with the exception of the SB counts 
below RUO' 17 erg cm' 2 s' 1 , where z > 0.6 star-forming 

17 See http://www.cfa.harvard.edu/hea/genx/ and 
http://hea-www.cfa.harvard.edu/SMARTX/ for more information about 
Generation-X and SMART-X, respectively. 


galaxies dominate. These galaxies are expected to produce 
X-ray emission from high- and low-mass X-ray binaries, 
hot gas, and young stars. By contrast, passive early-type 
galaxies, with X-ray emission produced by LMXBs and 
hot gas, make up only small fractions (»5%-20%) of the 
SB and HB number counts near the survey flux limits. These 
trends are well described by models that assume that (1) the 
X-ray power output for late-type galaxies scales with SFR 
out to z 2 with no evolution in the scaling relation and 
(2) the X-ray power output in early-type galaxies is due to a 
non-evolving hot-gas-emitting component and an evolving 
(by (1 + z) 2 ) LMXB component that both scale with stellar 
mass. Therefore, the rapidly rising normal-galaxy number 
counts can be attributed primarily to the evolution of the 
physical properties of galaxies (i.e., SFR and stellar mass) 
and not X-ray scaling relations. 

7. Extrapolation of our number counts to lower fluxes suggests 
that an X-ray observation reaching flux levels just below 
those probed by the ~4 Ms CDF-S (at SB fluxes of 
R*(3-5) x 10' 18 erg cm' 2 s' 1 ) would result in an X-ray sky 
with normal galaxies dominating the source density at the 
sensitivity limit. Since the normal-galaxy number counts 
are now comparable with AGNs at the R»4 Ms flux limit 
and are expected to continue to rise rapidly at fainter fluxes, 
more sensitive surveys will have relatively large yields in 
increased source densities. We show that a 10 Ms CDF-S 
would yield 020- 1080 total sources and would result in 

a limiting SB source density of ~46,000 deg' 2 with normal 
galaxies dominating the number counts (providing R»57% 
of the total number counts) at the flux limit. 
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APPENDIX 

RECOVERY FRACTION CORRECTIONS 

As noted in Section 2.2, the XI 1 CDF-S source catalog 
was generated following a two-step approach, which entailed 
(1) running wavdetect to form an initial list of candidate 
sources and (2) assessing the probability of detection using 
AE. Our methods described in Section 2 for computing number 
counts generally account for biases due to completeness and 
the Eddington bias provided that our catalogs are complete 
to a specific AE source selection probability. Therefore, we 
need to account for incompletenesses in the number counts that 
result as a consequence of our two-step cataloging approach. 
In this section, we summarize our approach for computing this 
correction, which has been implemented in Equation (5) as the 
quantity C' . 

We began by generating 200 mock CDF-S images in the four 
bandpasses (i.e., 800 images), which each contained 700 X-ray 


21 


The Astrophysical Journal, 752:46 (23pp), 2012 June 10 


Lehmer et al. 



Figure 17. Fraction of simulated sources satisfying the binomial probability selection criterion of Equation (1) that were also detected using wavdetect at a 
false-positive probability threshold of 10~ 5 (i.e.,/ rec0 v; filled circles with error bars) for the four bandpasses. Each color represents a different off-axis angle interval 
as annotated on the plot. Our best-fit models, as described in the Appendix, have been plotted as curves (see annotation); these models have been used to correct for 
incompleteness in our number-count computations in Section 2. 

(A color version of this figure is available in the online journal.) 


sources artificially implanted. These images were constructed 
as follows. Each source was assigned a random right ascension 
and declination within the area of the ~4 Ms CDF-S as defined 
by the exposure maps (see Section 3 of XI 1) and was given 
a random number of counts between ~ I and 100. The counts 
from each source were added to a blank image canvas using 
the marx ls (version 4.5) ray-tracing code. In this procedure, 
the PSF of each source was modeled assuming a roll angle and 
aim point that were selected from one of the 54 observations 
that make up the cumulative ~4 Ms CDF-S exposure (see 
Table 1 of XI 1). Selection of the roll angle and aim point for 
each source was done probabilistically with the probability of 
selection being directly proportional to the exposure time of 
each of the 54 CDF-S observations. Mock images for each 
of the four bandpasses were then created by adding background 
counts from the respective cumulative background images (see 
Section 7.1 of XI 1) to the image canvases. 

For each of the 700 sources in each of the 800 images we 
performed circular aperture photometry using a circular aperture 
with a radius encompassing ~90% of the EEF. The number of 
counts measured for each source was then compared with the 
number of counts needed for a source detection (see discussion 
of sensitivity maps in Section 2.1) to see if it would satisfy the 
criterion adopted in Equation (1) for source detection. We then 


18 See http://space.mit.edu/cxc/marx/ for marx simulator details. 


searched each of the 800 images using wavdetect at a false- 
positive probability threshold of 10~ 5 and made wavdetect 
source catalogs for each image. For each bandpass, we combined 
all 200 catalogs of 700 sources and computed the count rate 0 
and off-axis angle 0-dependent fraction of sources that satisfied 
our binomial probability selection criterion (i.e., Equation (1)) 
that were also detected by wavdetect (hereafter recovery 
fraction / recov = 1/C'). 

In Figure 17, we show the/ rec ov versus 0 for off-axis angle 
intervals going to 0 9 arcmin, where the image is roughly 

radially contiguous. We found that the data for / reC ov could be 
represented successfully using the analytic form / re cov(0, 0) = 
1/(1+ exp[-<5(0){0 - £(0)}]), where 8(6) and | (0) are fit- 
ting constants that vary with off-axis angle and bandpass. In 
Section 2.2, we made use of this relation in Equation (6), where 
we adopt C' = 1 //recov- We note that this correction has the 
advantage in that it is independent of the Bayesian priors dis- 
cussed in Section 2 and does not depend on the shape of the 
number-count distribution near and below the flux limits. 
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